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£H ■ Abstract 

d . 

We start with a short introduction to the roof concept. An elementary discussion 
' of phase-damping channels shows the role of anti-linear operators in representing 

their concurrences. A general expression for some concurrences is derived. We apply 
it to 1-qubit channels of length two, getting the induced foliations of the state space, 
the optimal decompositions, and the entropy of a state with respect to these 
channels. For amplitude-damping channels one obtains an expression for the Holevo 
qq , capacity allowing for easy numerical calculations. 
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^ : 1 Introduction 

p . 

The aim of this paper is in calculating some entanglement quantifying functions for a 
rather restricted class of quantum channels, in particular for all 1-qubit channels of 
length two. These channels are now well examined and classified due to the work of 
Fujiwara and Algoet pQ, King and Ruskai Ruskai et al [3J, Verstraete and Verschelde 
[3]. They include all extremal 1-qubit channels and some important doubly stochastic 
ones. An introduction to quantum channels is in Nielsen and Chuang, 



A channel, say T, is of rank two, if T(p) is of rank two for all density operators. In this 
case there are at most two eigenvalues of T(p) different from zero. It follows from trace 
preserving, that T(p) is characterized by its largest eigenvalue up to a general unitary 
transformation. Nevertheless, the treatment of a general rank two channel is beyond our 
present abilities if one asks for quantities like capacity or concurrence. 

A completely positive map is of length two, if it can be written down with two Kraus 
operators, but not with one. 

Let us now repeat some definitions and properties of the quantities to be treated. An 
ensemble is a finite number of density operators, every one given with a definite 
probability or weight, 

8 = {pi, . . . ,p m ;pi, ■ ■ ■ ,p m }, (1) 

such that the sum of the positive numbers pj is one. The density operators are states of 
a physical system. We may think of a quantum alphabet with quantum letters pj, and of 
a quantum message which chooses randomly the quantum letter pj with probability pj. 

The weighted sum of (^Q) is the average density of the ensemble, 
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While the quantum letters of £ are states of a physical system, the average (J2J is a 
"non-commutative probability measure" . Its role is quite similar to a "classical" discrete 
probability distribution. The extractable information per quantum letter of the quantum 
message is bounded by Holevo's jS] 

X{£) ■= S(p) - ^2pjS(pj) = ^PjS{p || pj), p = av[£] (3) 

Here S(.) is the von Neumann entropy and S(. || .) the relative entropy. 
Remark 1: The important inequality 

X (T(£)) < X (£) 

follows from the monotonicity of relative entropy which can be proved for trace 
preserving, at least 2-positive maps. Counter examples for just positive maps seems to 
be unknown. There are two recent reviews on the finer properties of relative entropy, one 
by Petz going beyond complete positivity, and one by Ruskai jSJ. 

Let T be a positive trace preserving map acting on the density operators of 7i. The 
application of p — > T{p) generates, letter by letter, a new quantum message belonging to 
the new ensemble 

T(£) = {T(p 1 ),...,T(p m y,p 1 ,...,p m } (4) 
The 1-shot or Holevo capacity of a channel T is the number 

C«(T):=inaxx(T(£:)) (5) 

We are aiming at a slightly more delicate expression. To see its significance we try to 
perform (J5J) in two steps. At first we let £ run only through the ensembles with a given 
average density p, postponing the search for the maximum, i. e. the Holevo capacity, to a 
later time. In doing so we first define 

H(T;p):= max X (T(£)) (6) 

a,v[t\=p 

so that the maximum is to compute with respect to all ensembles having p as its average 
density. We get a function, depending on p (and T) only. If T is a partial trace onto a 
sub-algebra, © has been called entropy of p with respect to the sub-algebra by Connes et 
al see also Narnhofer and Thirring ^Uj. In this sense © could be called entropy of p 
with respect to T. 

Schumacher and Westmoreland JJ| relate © to the efficiency of quantum channels. 
They denote the quantity in question by x*(p) an d they call an ensemble saturating © 
an optimal signal ensemble. 

It is important for our purposes to rewrite (jHJ) by the help of (j3J) as 

H(T;p) = S(T(p))-E(T;p) (7) 

where E is an entanglement measure given by 

E(T;p)= min 5>S(T(p,)) (8) 
2^PkPk=p 

and the minimum is to compute over all convex decompositions of p. If T is a partial 
trace in a bipartite quantum system, (jHJ) is the entanglement of formation introduced in 



2 



Bennett et al followed by the remarkable papers of Wootters, ^3]) and of Terhal 
and Vollbrecht |14j . Other examples are in Benatti et al |15j . |16) . In ^7] some relations 
between the quantity © for quite different channel maps are pointed out by the same 
authors. 

It is well known that E(T ; p) is convex in p. Indeed, E is written as the convex hull of 
the function S(T(p)) in ©. The concavity of the entropy allows to perform (jHJ) over the 
pure decompositions of p without changing the outcome of the minimization. Thus 

E(T;p) =mmJ2PjS(T(TTj)), p = J2Pj^j (9) 

where all the irk are pure and the minimum is to perform with respect to all pure 
decompositions of p. An ensemble of pure states saturating © will be called an optimal 
ensemble. 

One should also notice: The equality between the two variational problems (|HJ) and © is 
due to the concavity of S(T(p)). If we perform similar computations, however with a 
function not being concave, the two problems are essentially different. For certain 
channels, including all rank two and length two 1-qubit channels, the minimization @ 
will be solved in the present paper. 

The next section is a short introduction to the roof concept. We shall discuss how to 
estimate © from below by convex functions and from above by roofs. 

Next we demonstrate the procedure for the phase-damping channels, where most things 
are now well understood. It follows the computation of concurrences and of E for some 
rank two channels. We give quite explicit computations to see the dependence of 
concurrences, foliations, and optimal ensembles from the Kraus operators. 

A quite interesting observation is the following: Let us call Kraus module of T the linear 
span of the Kraus operators of T. Choi ^S] has shown that irreducibility of trace 
preserving cp-maps is a sole property of the Kraus module. Here we prove for length two 
and rank two channel maps that their optimal ensembles and, hence, their foliations are 
equal if they belong to the same Kraus module. It is not probable that this remarkable 
feature will survive for more complex channels. Nevertheless it seems worthwhile to ask 
for similarities of channels with identical Kraus modules. 

It is a further consequence from our calculations that the foliations deform continuously 
by changing the entries of the Kraus operators for the class of channels considered. The 
foliations are more coarse properties of rank two cp-maps than concurrences. 



2 Roofs 

There are some general features in the optimization problems we are interested on. They 
constitute the ground floor for more refined investigations. 

Let us abstract from the specific values given at the pure states in © and let us start 
with an arbitrary real valued and continuous function g(ir) on the set of pure states. 
Assume, we like to extend g to a function defined for all density matrices. Clearly, there 
are many and quite arbitrary solutions for the problem. Let us denote by G one of these 
extension. 
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To place a first restriction, we require the extension to be "as linear as possible" . The 
requirement can be made precise as following. Let p be a density operator. If there exist 
a pure decomposition of p such that 

G (p) = Y,Pi9(Kj), P = ^>2Pjnj, (10) 

we call the decomposition p-optimal or simply optimal for G. We also call a pure 
ensemble 

£ = {it i, . . . ,n m ;pi, . . . ,p m }, Pj > 0, 

G-optimal if 

A function, G, which allows an optimal decomposition Q1U|) for every state p, I call a roof 
or, more literally, a roof extension of g. Roof extensions reflect the convex structure of 
the state space and they are, in a well defined way, "as linear as possible" . It is not easy 
to gain good examples of roof extensions in higher dimensions. In general, however, 
there are a lot of them for a given function g. 

Let us now consider two further possibilities to extend g from the pure ones to all states: 
We may require the extension to be either concave or convex. In this spirit we call a 
function on the state space a convex extension of g if the extension is a convex function. 
Similar we speak of a concave extension of g if the extension is concave. 

Lemma 1. Given a convex, a concave, and a roof extensions of g. Then 

(P) (11) 

for all density operators p. 

The proof is almost trivial: One chooses a p-optimal decomposition (|l(Jj) . Then the very 
definition of a convex (concave) function establishes (jll|) . 

It follows that, given g, there can be only one convex roof extension, "the" convex roof 
determined by g. The family of roof extensions of g has just one member in common 
with the family of its convex extensions. Similarly there is just one concave roof which 
extends g. 

As a matter of fact |19j . the convex roof with values g(ir) at the pure states is nothing 
than the solution of the variational problem which mimics HU • 

Gconvex,roof(/3) = m[n ^2 P j 9 j) , P = ^Pj^j ( 12 ) 

Here one has to run through all pure state decompositions of p. We get © by setting 
g(ir) = S(T(tt)) for pure ir. On the other hand, if we take the maximum in ()12|) instead 
of the minimum, we obtain the concave roof extension of g. 

Let us return for a moment to the more specific of evaluating ©. To calculate E(T p) 
amounts to construct the convex roof with the function ^(vr) = S(T(ir)). Looking at the 
roof property, there are some typical questions one should ask. For instance we may 
start by a set of parameterized mappings, T s , and we like to know whether they have, 
perhaps for some p, identical optimal decompositions. More literally, we ask for a pure 
decomposition of p which is optimal for every T s in an appropriate range of the 
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parameter s. In the next sections we like to convince the reader that this point of view is 
quite fruitful. To do so, a remarkable property of convex and concave roofs is to explain. 

Let G be a convex roof on the density operators of a Hilbert space of finite dimension d. 
At first we use convexity: Let us fix a density operator p. There is at least one 
Hermitian operator, say Y p , such that for all density operators ui 

G{p) = Tr Y p p, G(u) > Tr Y p uj (13) 

is valid. Now let us apply the roof property: There is a pure decomposition 

9 = ^Pj^j, Pj > 0, TTj pure, 

which is G-optimal. Thus 

$>Tr Y% = Tr Y p p = G(p) = ^PjGfa) 

Because of the inequality 13|) this can hold if and only if 

G(Tr j ) = TvY%, j = 1,2,... 

Let us now look at the convex set 

{ u : G(u) = Tr Y p uj } (14) 

By the very construction, G is convexly linear (affine) if restricted to this set. On the 
other hand, it contains p and, by the reasoning above, it contains every pure state which 
belongs to a G-optimal ensemble with average density p. Or, in other words, Q14|) 
contains all the pure states which can appear in an optimal decomposition of p. Let us 
collect all that in a theorem, jl^j - 

Theorem 2. Let g be a continuous function on the pure states. 

i) There exist exactly one convex roof extension, G, of g. 

ii) G can be represented by the optimization procedure \12\). 

Hi) There exist optimal pure decompositions for every density operator p. 

iv) Given p, G is convexly linear on the convex hull of all those pure states which can 

appear in an optimal decomposition of p. 

Let us call the convex hull of all pure states appearing in all possible optimal pure 
decompositions of p the optimal convex leaf of p. As proved above, and stated in the 
theorem, G must be convexly linear on every optimal convex leaf. 

Remark 2. The remark concerns property hi). The compact convex set of density 
operators enjoys a peculiarity: The set of their extremal points, i. e. the set of pure 
states, is a compact one. This allows to prove the existence of optimal decompositions 
IjlOJl by the continuity of g. Then, by a theorem due to Caratheodory, one deduces the 
existence of optimal decompositions with a length not exceeding d 2 , d the dimension of 
the Hilbert space which carries the density operators. It should be noticed that the 
compactness of the extremal points is an extra property. Counter examples are by no 
means exotic as seen by the set of trace preserving cp-maps. To get the conclusion hi) it 
suffices for g to be continuously extendable to the closure of the set of extremal states. 

Remark 3. If G is a convex roof and G is the sum of two convex functions G = G\ + G2, 
then G\ and G2 are convex roofs, and every G-optimal ensemble is optimal for G\ and 
G 2 . 
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We need one more definition: We call a convex roof, G, flat if it allows for every p an 
optimal pure decomposition 

P = SPi^ g (p) = Y.Pi G ^i) 

such that 

G(-K 1 ) = G(TT 2 ) = ... = G(-K j ) = ... 

If this takes place, every p is contained in a convex subset, generated by pure states, on 
which the roof is not only convexly linear but even constant. The merit of flat roof, say 
G, is in the nice property that every function of G, say /(G), is again a roof. 

As a matter or fact, the convex roofs we are considering in the following enjoy even a 
stronger property: They are constant on the convex leaf of every p, i. e. g is constant on 
the pure states of every pure ensemble which is optimal for G. 

3 Phase damping channels 

Let us consider some particularly simple examples of 1-qubit channels, the family of 
phase-damping channels, to see what is going on in applying concurrences, ^F3\, to 
1-qubit channels according to (20]. For the most symmetric channel of the family, with 
z = in l|15|l. the theorem below and the insight into the foliation of the state space are 
due to Levitin |2T] . 

Let \z\ < 1 be a complex number. Define the map T z by 

x= (xoo xoi\ ^ Tz{x) = ( y zx 01 \ 

\Xio X\l J \z*x w XXI J 

The application of such a map does not change the pure states |0)(0| and |1)(1|, and 
there are no other trace preserving and completely positive 1-qubit maps with this 
property than those given by ((T5|) . 

Before starting the calculation let us have a look on a bundle of parallel lines which 
foliates the state space. Geometrically, the 1-qubit state space can be represented by the 
Bloch ball. The Bloch ball is the unit ball sitting in the "Bloch space", that is in the real 
Euclidean 3-space of all Hermitian matrices of trace one. The "Pauli coordinates", Xj, of 
a matrix are read off from 

X = - (x l + XX&X + X2&2 + 

With xo = 1 the real coordinates x\, X2, X3, parameterize the Bloch space, and in this 
context they are referred to as "Bloch coordinates" . Finally, the Bloch ball is the unit 
ball with respect to the Bloch coordinates. 

Let us return to the phase-damping channels. The line x\ = X2 = 0, i. e. the X3-axis, 
remains point- wise fixed under the mappings (|15|) . On the intersection of the line with 
the Bloch ball E(T Z ; p) must be zero. One aim is to show that E(T Z ; p) is constant on 
the intersection of the Bloch ball with every line which is parallel to the line x\ = X2 = 0. 
These lines are given by fixing the values of x\ and X2 and letting X3 free. Equivalently, 
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such a line can be given by fixing xoi in Ijl5|) . Now we can write down a very simple 
convex roof, the restriction of the function 

X i y Jx\ + x\ = 2 |xoi| 

onto the state space. Indeed, this function is convex on the Bloch space. It is even a 
semi-norm there. And it is, trivially, a flat roof: It is constant on every line parallel to 
the X3-axis of the Bloch space. The intersection of a given line with the Bloch ball is 
either empty, or touching the ball in just one point, or is a line segment with two pure 
states as end points. In the latter case, the pure states of the segment are 

V xio p J \xio l-pj 

with 2|xoi| < 1. For 2|xoi| = 1 the line touches the Bloch ball at one pure state. 

With a flat roof one can build other roofs just by taking a function of it: If / is a real 
function on the unit interval, /(|xoi|) is again a flat roof. By an appropriate choice of / 
we shall find the form of E(T Z ; p). 

To do so we need to compute the determinant 

detT 2 (X) = xoo^ii — zz*xq\X\q = detX + (1 — zz*)xq\Xiq 

Taking X pure, say Xjk = OjO^, only the second term is different from zero and we 
remain with 

detT z (X) = (1 - |z| 2 )|a ai| 2 , x jk = aja* k (17) 
Using ideas from ^3] an d |H] we define the concurrence of T z for Hermitian X by 

C(T Z -X) := v /(l-| z |2)(xf + xi) (18) 

The definition differs from the one used in [20] by a factor two. It influences, here and 
later on, the appearance of some equations. The concurrence is a semi-norm on the 
Bloch ball, and, if restricted to the Bloch ball, it is the unique convex roof satisfying 

C(T Z ; tt) = 2y / detT,(vr), tt pure (19) 
Following again ^2] and ^3] we introduce 

h(x) = —x log x — (1 — x) log(l — x) 
and, using ad hoc notations, 

h 1 (x) = h(^) 

and 

h 2 (x) = h 1 (y), y = \/l-x 2 

Theorem 3. For all \z\ < 1 and all density operators p 

E(T z ;p) = h 2 (C(T z ;p)) (20) 

holds. It is E{p\) = E(p 2 ) for two density operators, p\ and p 2 , if and only if they have 
equal distances to the x^-axis of the Bloch space. The pairs of optimal pure states are 
given by U6\) . 
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Proof: We already know that (j^Ojl is a roof with the desired values at the pure states. We 
only have to show that it is convex. Then, by the uniqueness theorem, we are done. One 
calculates the first and the second derivative of h%, assuming log = In. At first we get 

, x x l + y 

hn(x) = — In 

2K ' 2y 1-y 

For x > and y > we find h' 2 > 0. One further obtains 

h>i( X )= l - + y -+ y - + y - + ... 
2y ' 3 5 7 9 

and h'2 > proves the convexity of h 2 ■ Let C (p) be a convex function on the state space 
( - or on another convex set - ) with values between and 1. Let us denote by a dot the 
differentiation of C in an arbitrary direction. Then 

h 2 {Cy = h' 2 \C)CC + h' 2 {C)C 

By the convexity of C one gets C > 0, and we have seen h! > on the unit interval. 
Thus, the second term is not negative. As we know h" > 0, we have shown the convexity 
of the function (|2U[). and we done. 

Lemma 4. Let C be a convex function with values in the unit interval, defined on a 
finite dimensional convex set. Then h%{C) is convex. 

It is remarkable that the whole set of phase damping channels induces a single foliation 
of the state space: The foliation is a property of the Kraus module belonging to the 
channels. The single foliation forces the concurrences to differ by a factor only if z is 
changing. 

Now we have to add a further structural element as a guide in treating more general 
1-qubit channels. What we have in mind is, up to a contracting factor, a reflection on 
the plane 23 = 0. In Bloch space a reflection is not a proper rotation. Its functional 
determinant must be negative, enforcing its implementation by an anti-linear operator in 
the 1-qubit space. The anti-linearity is unavoidable. 

Let us define an anti-linear operator -d z by 



'■(Z) = (21> 

The operator is an Hermitian one, i. e. it satisfies 

ifaWzlfc) = (<h\#z\<h) 

for all pairs of vectors. Being anti-linear, 1? acts to the right, but not to the left. With 
two arbitrary vectors, 

">-(£)■ ">-(£)■ 

we get the following relation. 

(a\# z \b) = VI - zz* {aoh + ai&o)* (22) 
In combination with (|17|) we obtain the equation 

4detT z (|a)(a|) = \{a\$ z \a)\ 2 (23) 
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which, together with (jl9fl . can be written 

C(T z ;\a)(a\) = \(a\& z \a)\ ( 24 ) 

Because the concurrence is a convex roof, we can now return to one of the properties of 
such function: 

C(T z ;p)=mm^2p j \(^Wz\^ j }\ (25) 

where the minimum runs through all possible ways of representing p as a convex 
combination 

P = £Pil#j)(&l 

The minimum will be attained by choosing the decomposition with a pair (jlfij) of pure 
states. With the same optimal decomposition of p the optimization problem for E(T Z ; p) 
can be saturated. (The diagonal entries of p must 1/2.) 

Remark 4- Let us assume we like to solve for the phase-damping channels the 
variational problems (|SJ), but we like to replace the minimization by the maximization, 
resulting in concave roofs. Now the foliation for the maximization is given by the 
intersection of the planes perpendicular to the X3-axis with the Bloch ball. The foliation 
is the same for all z. The intersection of a plane with a ball contains a whole circle of 
extremal states. Hence there are very many different optimal decompositions for given 
mixed p. 



4 Concurrence 

Let us now discuss some generalities for rank two channels and let us see, where the 
difficulties are. A quantum channel, T, is of rank k if the maximal rank of all pictures, 
T(p), is k. 

Let T be of rank two. We assume in addition that T maps into a 1-qubit space. Then 
the determinant of T(p) is the product of the eigenvalues of T{p). We define the 
concurrence, C(T ; p), of T to be the convex roof which attains at pure states the values 

C(T;vr) = 2^/det T z (vr), vrpure (26) 

completely similar to (jl9|) . Being convex with values between and 1, we can literally 
repeat the construction as in Q2U|) . The result is a convex function which coincides at 
pure states with E. But in general we do not know, wether (j26|) is a flat roof. Hence, we 
only can conclude 

E(T;p)>h 2 (C(T;p)) (27) 

If, however, (|26|) is a flat roof, then the right hand side of i|26j) is a flat convex roof and 
we have two convex roofs agreeing on pure states. Then equality must hold. 

Lemma 5. If C as given by \2fy) is a convex roof, then we have 

C flat E(T;p) = h 2 {C(T-p)) (28) 

It is useful to know whether the concurrence of a channel map T is a flat roof. The 
following theorem gives a whole class of them. 
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Theorem 6. Let i? be an anti-linear hermitian operator and define C$ as the convex 
roof extension of 

C,(|$M) = |(0W>l (29) 
for all pure states ir. Then the convex roof C$ is flat. It is 

C#(p) = max{0,A 1 - £>i } (30) 

where \\ > A2 > . . . are the ordered eigenvalues of 

{pVHp&p 1 ' 2 ) 1 ' 2 

Corollary 7. If there exists an anti-linear hermitian i9 such that 

4detT(7r) = Tr 7n?7n?, 7rpure (31) 
is valid for all pure density operators then 

a a (p)=C{T;p) (32) 

is a flat roof and (|28|) is valid. 

This theorem is proved in [22] . It provides a certain application of the methods of 
Wootters and others, see Wootters |23| . As already mentioned, with T the partial trace, 
E(T ; p) is the entanglement of formation. In the 2-qubit system i9 is the Hill- Wootters 
conjugation. 

For 1-qubit maps ()30j) reads Ai — A2 and one can become more explicit. Abbreviate 

t r 1/2 1/2x1/2 

i = (pi P2P1 ) 1 

Taking the trace of the characteristic equation of £ results in 

Tr£ 2 + 2det£ = (Tr£) 2 
and the squared sum of the eigenvalues of £ becomes 

(Ai + A 2 ) 2 = Tr(p lP2 )+2det£ 

Combined with the relation 

(Ai-A 2 ) 2 = (Ai + A 2 ) 2 -4det£ 

it yields 

(A x - A 2 ) 2 = Tr ( PlP2 ) - 2 v /(det / oi)(det / 9 2 ) 
Substituting p\ = p and p 2 = ^P"& provides, [2*2*] . 

G & (p) 2 = Tr (p-dpd) - 2(det p) det v 7 ^ (33) 
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5 1-qubit channels of length two 



In this section we like to show that the corollary to the preceding theorem applies to 
1-qubit channels of length two. As we shall see, the existence of an anti-linear hermitian 
9 fulfilling ()31|) does not depend on trace preserving. 

Let A and B be two linear independent operators on a 2-dimensional Hilbert space and 

^ _ ( a oo a oi \ B = ( ^ 01 
\a 10 an) 1 \ho hi 

their matrix representations with respect to a reference basis. 

T(X) = AXA^ + BXB^ (34) 

is a completely positive map of length two. The following has been proved in [2U] by 
straight forward computation. 

Theorem 8. There is an anti-linear hermitian operator $a,b such that for all pure 
states 

4detT(|o)(a|) = \(a\$A,B\b)\ 2 

where T is given by \34\) - 

If such an operator exists, it is determined by T up to a phase factor only. The 
ambiguity is a natural one due to a geometric phase. 

To describe one can introduce its matrix representation 

tfA£( a °w aoo i +aoia n (35) 

' \aij \a 10 a^ + a n alJ 
and express the matrix elements by those of A and B. 

«oo = 2(&ioa o - aioboo)*, = 2(a i&n - 6 ian)*, (36) 

«oi = "io = (aoofrn - an&oo + aoi&io - aio^oi)* (37) 

Up to a factor, due to another normalization of the concurrence, this is agreement with 
|2Uj . Clearly, $ must be skew symmetric in the matrix entries of the Krause operators, 

$A,B + $B,A = 0. 

Comparing the equations above with phase-damping channels one gets $a,b = —^z- 

One may expect a more transparent representation than (|36j) and (|37|) . This is possible 
by the spin-flip operator Of, the "fermion conjugation" for one qubit, 



After calculating 



A^QtB 



(frioaoo — aio^oo)* (aoo^ii — aio^oi)" 



/ V (ooibio ~ a n fo oo)* (ooi&n - aii&oi)" 
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and comparing that expression with and (|37|) . we get 

$ AB = A^9 f B - B ] 9 f A (38) 

Remember 9j = —Of = 9J 1 to see that (j38[l is an anti-linear hermitian operator - as it 
should be. Now assume a transformation of the Kraus operators according to 

A, B — ► CiAC 2 , C X BC 2 (39) 

Because of 

C\e f Ci = (detCi)*0/ 
the anti-linear operator must change as follows: 

#a,b — > (detd)* cl# AjB C 2 (40) 

That it is of value to classify super-operators up to a transformation (|39|) is also seen by 
the results of Verstraete and Verschelde, (I]. 

Let T' be another length two cp-map with Kraus coefficients A' and B' , and let us 
assume a linear dependence 

A' = nooA + noiB, B' = {i 10 A + fi u B (41) 

By the help of (|38[). or by observing that (|3f))) and (|3*7|) can be expressed by determinants 
in the coefficients of the Kraus operators, one can reproduce the relation 

"&A',B' = (W)0/*ill - Moi^io)*^a,b (42) 

Why is this interesting? It shows that our procedure associates, up to a scalar factor, to 
every pair of operators, chosen from the linear span of A and B, the same anti-linear 
operator. With other words, to every 2-dimensional Kraus module a 1-dimensional linear 
space of anti-linear hermitian operators is attached. 

Remark 5. Regard the Kraus modules for the 1-qubit channels as the points of the 
second Grassmann manifold of the space of linear operators. By attaching the multiples 
of ??a,_b to the corresponding points one gets the line bundle. It is dual to the 
determinant bundle as one can deduce from what follows. For the time being, we shall 
not follow further this way. 

An observation, related to (|12*|). is 

A' (8) B' - B' (8) A' = (/xoom - M01M10) (A® B — B® A) (43) 

If we apply the operator identity (|43j) to a \aa), we get an anti-symmetric 2-qubit vector. 
There is, up to a factor, only one such vector and the yet unknown factor must 
transform as in ()43l) . Performing the calculations one gets 

(A ® B - B <g> A) \aa) = ^ (a\ti A , B \a)* (|01) - |10» (44) 

Remember that for a single channel, T, only the absolute value of the expectation value 
at the right hand side is relevant. 
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Two channels may be called "unitary equivalent" if 

T'(p) = U 1 T(U 2 1 pU 2 )U{ 1 
with two unitaries XJ\ and U 2 ■ For trace preserving channels we can assume 

A =(T :,)• B =(L "?) w 

up to unitary equivalence, |Hj . Then $ becomes diagonal, 

& A , B ( ao )=2( (Wor \ (46) 

Vai/ V(-ooionai) */ 

Abbreviating 

yo = 26i aoo, 2/1 = 26 i«n 

we can write 

PooVoVo -pioyoy* 



V -poiViVo PooViy'i 
It follows 

Tr ptfptf = Poo2/o2/S - Poi2/oyi - Pio^o + Poorly* 

On the other hand, 

2 det p det VlP = 2 \y y 1 | (pooPn - P01P10) 
Pasting all things together, we get from 

C 2 = (pool 2/0 1 - P11M) 2 + 2\y yi\p Q ipio - Poi2/o2/i - Pi 2/o2/i 

We now choose the square roots of yoy% and y$yi such that their product becomes 
positive. Then we can write the remaining terms above as follows: 

-{pm\Jyoy\ - pm\Jylyi) 2 

so that, finally, we see: 

C(Tp) 2 = AL 1 (p) 2 + 4L 2 (p) 2 (47) 
were L\ and L 2 are real valued and linear in the entries of p, 

Li(p) = Poo\bi a o\ - pii|&oiaii| (48) 



L 2 (p) = i(poi^Jhoaooboi a *i ~ Pio\/b* a^ o b 01 au ) (49) 
and we have to choose the signs of the roots according to 

^/&ioaoo&oi°ii \f b lo a oo b oiaii > (50) 

The result compares well with the more symmetrical case of the phase-damping channels: 
C is the square root of a positive semi-definite quadratical form. Geometrically, the 
points of constant concurrence and, hence, of E are ellipse-based cylinders. 

In the non-degenerate case none of the two linear forms vanish identically. The foliation 
of the state space at which the concurrence and also E{T ; p) remain constant are given 
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by the straight lines which are the intersection of the two families of planes L\ = 
constant, L2 = constant', in the Bloch space. There is just one straight line at which the 
concurrence is zero. It goes through the two pure states which are mapped onto pure 
states by T. Up to normalization these two pure states belong to the vectors 



a 



^ , frioaooao — &oi a nai 
They are mapped by T to vector states of the form 



Let us now shortly look at the degenerate case in which 601610 is zero in (|45j) ■ The 
amplitude-damping channels are well known examples. They can be defined by the 
Kraus operators 

with < p < 1. The action of T is described by 

Poo Pol \ , ( Poo + (1 -p)pu V 1 - PP01 



(|47j) becomes 



Pw P11J \ V 1 ~ PP10 PPn 



C(T;p)=2 y J P (l-p)p 11 (52) 

The two families of planes degenerate to one family, the planes perpendicular to the 
3-axis of the Bloch space. The foliation dictates the behavior of C(T ; p) and E(T ; p). 

One observes that, given pn, the entropy S(T(p)) becomes maximal if the off diagonal 
entries of p vanish. Therefore, if p' is the diagonal part of p, We get 

E(T;p)=E(T;p'), H(T ; p) < H(T ; p) 

and to obtain the Holevo capacity it suffices to consider diagonal density operators only: 

C^ 1 ' = max H(T ; p'), p diagonal 

Writing r for pn, such that < r < 1, we can rewrite the capacity as follows: 

C«(D = max \h(pr) - »( 1 = ^ ~ f S E **" )] (53) 

0<r<l z 

Because H(T ; p) is a concave function, (JSHJ is a concave function of r and, obviously, not 
degenerate. Therefore, for any given < p < 1, there is exactly one value ro of r at 
which the maximum in is attained. 
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